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Dependence of the optical brightness on the gamma and X-ray 
properties of GRBs 
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The Swift satellite made a real break through with measuring simultaneously the gamma X-ray and optical 
data of GRBs, effectively. Although, the satellite measures the gamma, X-ray and optical properties almost in 
the same time a significant fractions of GRBs remain undetected in the optical domain. In a large number of 
cases only an upper bound is obtained. Survival analysis is a tool for studying samples where a part of the cases 
has only an upper (lower) limit. The obtained survival function may depend on some other variables. The Cox 
regression is a way to study these dependencies. We studied the dependence of the optical brightness (obtained 
by the UVOT) on the gamma and X-ray properties, measured by the BAT and XRT on board of the Swift 
satellite. We showed that the gamma peak flux has the greatest impact on the afterglow's optical brightness 
while the gamma photon index and the X-ray flux do not. This effect probably originates in the energetics of 
the jet launched from the central engine of the GRB which triggers the afterglow. 



1. INTRODUCTION 

A significant achievement of the Swift satellite is 
the simultaneous detection of the physical properties 
of the gamma ray bursts in the gamma, X-ray and op- 
tical domain, measured by the BAT, XRT and UVOT 
instruments on board of the satellite. 

Following the alert given by BAT the satellite starts 
to slew and after reaching the position of the burst the 
XRT and UVOT make measurements in the X-ray and 
optical domain, respectively. 

Although, a significant fraction of the bursts is de- 
tected by the XRT as well, it is not the case with the 
UVOT where at a remarkable fraction of the events 
only an upper limit of the optical brightness is ob- 
tained. 

From theoretical point of view the measured opti- 
cal and gamma properties may be given by completely 
different phenomena, their observational relationship, 
if there is any, would be an important constraint for 
the possible models. To study this relationship it 
would be a serious bias if we take into account only 
those cases where all properties, i.e. gamma, X-ray 
and optical, are measured. 

Survival analysis is a way to make use the informa- 
tion which is inherent in the value of the upper bound 
of the optical brightness. Cox regression is a tool for 
studying the dependence of the survival function (a 
result of the analysis) on some background variables, 
the covariates (gamma and X-ray properties in our 
case). 

In the following we use Cox regression to study 
the dependence of the distribution of the UVOT de- 



tected optical brightness on measured gamma and X- 
ray properties. 



2. MATHEMATICAL SUMMARY 

Let we have a t stochastic variable with f(t) prob- 
ability density. The S(t) survival function is defined 
by 



L 

J f(t')dt' = F(t) = l-S(t) 



(1) 



where F(t) means the probability distribution func- 
tion. Actually, the S(t) s urvival function is its com- 
plement (F(t) + S(t) = l). iKaplan and Meyer I [l958j 
showed that S(t) can be estimated bias free even in the 
case when some of the values in the t\, ■ ■ ■ ,t n ob- 
served sample are only lower bounds (censored). The 
ratio of f(t) to S(t) is called the hazard function: 

f(t) S'(t) d 

m = W) = = ~dt l09[s{t)] (2) 

The h(t) hazard function characterizes the risk that 
in the [t, oo] range (S(t) gives its probability) an event 
will happen in the [t,t + dt] interval (its unconditional 
probability is f(t)dt). The hazard function may de- 
pend on bac kground var iables (the covariates). The 
Cox model (jCox I |l972j ]) assumes that this depen- 
dency can be written in the form of 



log[h(t)] = a(t) + Bxxx + B 2 x 2 



(3) 
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Figure 1: Distribution of the GRB positions according to 
the Galactic latitude. Note the depression in the 
distribution close to low latitudes due to the foreground 
extinction of the Galactic dust. Shadowed region marks 
the excluded area. 



where xi, x%, . . . , x m are the covariates while the a(t) 
arbitrary function and the B±, B2, ■ ■ ■ , B m constants 
have to be determined during the procedure of the Cox 
regression. If all these constants are equal to zero the 
a(t) function is identical with the logarithmic hazard 
function. The value of the constants characterize the 
strengths of the influence of covariates on the hazard 
and, consequently on the survival function. 



Table I Spearmann's correlations (numbers in bold face 
mark significant correlations) 





T90 


Flu 


Peak 


Pind 


Xflu 


V 


T90 corr. 
sign. 


1.000 


.646 

<.001 


-.066 
.149 


.155 

.001 


.408 

<.001 


-.024 
.646 


Flu corr. 
sign. 


.646 

<.001 


1.000 


.539 

<.001 


-.148 

.001 


.471 

<.001 


-.204 

.000 


Peak corr. 
sign. 


-.066 
.149 


.539 

<.001 


1.000 


-.282 

<.001 


.141 

.018 


-.289 

<.001 


Pind corr. 
sign. 


.155 

.001 


-.148 

.001 


-.282 

<.001 


1.000 


.049 
.406 


.004 
.944 


Xflu corr. 
sign. 


.408 

<.001 


.471 

<.001 


.141 

.018 


.049 
.406 


1.000 


-.109 
.080 


V corr. 
sign. 


-.024 
.646 


-.204 

<.001 


-.289 

<.001 


.004 
.944 


-.109 
.080 


1.000 



Table II Estimated B coefficients in Eq. Q (number in 
bold face means that the corresponding B differs from 
zero significantly.) 





B 


Wald 


df 


sign. 


logT90 


.732 


6.346 


1 


.012 


logFlu 


-.809 


4.099 


1 


.043 


logPeak 


1.886 


25.438 


1 


.000 


Pind 


.058 


.055 


1 


.815 


logXflu 


-.003 


.001 


1 


.973 



3. DESCRIPTION OF THE DATA 



We used for the present analysis the data avail- 
able in the Swift table 1 recorded until the date of 
03/03/2012, in particular, the V magnitude as a de- 
pendent variable, Duration, Fluence, Peak flux, Pho- 
ton index and early X-ray flux, as covariates in the 
analysis. Except of the Photon index we used loga- 
rithmic values in order to suppress the impact of the 
outliers on the regression. 

Since the optical brightness of the GRB afterglow 
is seriously dimmed by the foreground Galactic ex- 
tinction we excluded the cases with the latitude of 
|6| < 15°(|sin6| < 0.26). 

In the case if no afterglow was observed a lower 
bound in the stellar magnitude (upper bound for the 
observed brightness) was obtained. One can infer 
from Fig. Q] that at low latitudes the depression is not 
present in the distribution of the cases where only a 
lower V magnitue bound was determined. 



1 http://swift.gsfc. nasa.gov. /docs/swift /archive/grb.table 
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4. COX REGRESSION 



Before we make the Cox regression running we have 
computed the bivariate correlations between the vari- 
ables included in the analysis. In this procedure we 
have taken into account only those cases where both 
variable used for computing the correlation had mea- 
sured values. To overcome the problem with the out- 
liers we computed the Spearman's rank correlation 
which is not sensitive to these. 

In Table U we used all cases having measured val- 
ues in both variables used in the analysis, pairwise. 
We marked with bold face where the correlation co- 
efficients differ significantly from zero. As we are ap- 
proaching the detection limit of UVOT, however, only 
a lower magnitude limit is obtained in an increasing 
number of cases (see Fig. [2]). Of course, cases having 
only an lower limit in V have to be excluded from com- 
puting the correlation of this variable with the other 
ones. 

Results of the Cox regression are summarized in 
Table [TT1 Seemingly, the duration, gamma fluence and 
peak flux has a significant impact on the distribution 
of the optical brightness, while the gamma photon in- 
dex and the early X-ray flux do not. 
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Figure 2: Distribution of the afterglows' brightness in 
visual magnitude. Note the increasing fraction of 
censored data as we are approaching the detection limit. 



less, it is not the case with the gamma photon index 
and the early X-ray flux. 

The reason for the impact of some gamma proper- 
ties on the optical brightness is probably lying in the 
energetics of the jet launched from the central engine 
of the GRB which triggers the afterglow in the sur- 
rounding interstellar matter. 
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5. CONCLUSIONS 



We performed Cox regression in order to look for 
the impact of the gamma and X-ray properties of the 
GRBs on the afterglows' optical brightness. This ap- 
proach is necessary since in a significant fraction of 
cases only an upper bound of the optical brightness 
(lower bound in the V magnitude) can be determined. 

The analysis demonstrated that among the B coef- 
ficients in Eq. ([3]) belonging to the duration, fluence 
and peak flux differ significantly from zero. Neverthe- 
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